AD-A277  902 


NUWC41PT  Technical  Report  10,405 
1  October  1993 


A  Predictive  Parametric  Model  for  Volume 
Scattering  Strength  in  the  Deep  Ocean  With  an 
Application  to  Nighttime  Near-Surface  Scattering 
in  the  North  Sargasso  Sea  Below  16  kHz 

R.  Alfred  Saenger 
Wendy  L  Wachter 

Environmental  and  Tactical  Support  Systems  Department 


OTIC 

ELECTE 
APR  0  7  1994 


E 


94- 10486 


Naval  Undersea  Warfare  Center  Detachment 

New  London,  Connecticut 


Approved  for  public  release;  distribution  is  unlimited. 


94  4  6  013 


PREFACE 


This  document  was  prepared  under  Project  No.  D20Q27,  'T^uld-Static 
Stmar  Program  (MSS),"  Prinapal  Investigator  R.  Quisdan  (Code  31 12)  and 
Program  Manager  C  Mascm  (Code  304).  The  Sponsoring /^vity  is  thie  Rnogram 
Executive  C)ffice*Undersea  Warfare  (FEO-USW),  Project  Manager  Darrell  Spires 
(Advanced  Systems  and  Technology  OfSce-Cl)  under  Program  Hement  63SS3N. 

The  technical  reviewer  for  this  report  was  R.  CSmsdan  (Code  3112). 


Reviewed  and  Approved:  1  October  1993 


B.  F.  Cole 

Head,  Environmental  and  Tactical  Support 
Systems  Department 


REPORT  DOCUMENTATION  PAGE 


FomAppiovtd 
Om  No.  0^4-0188 


1.  AGENCY  USE  ONLY  (Lmv*  Blank)  2.  REPORTDATE 

1  October  1993 


4.  TTTLE  AND  SUBTITLE 

A  Predictive  Parametric  Model  tor  Volume  Scattering  Strength  to  the  Deep 
Ocean  with  an  Application  to  Nighttime  Near-^rface  Scattering  in  the 
North  Sargasso  Sea  Below  16  kHz 


6.  AUT>OR(S) 

R.  Alfred  Saenger  and  Wendy  L.  Wachter 


3.  REPORT  TYPE  AND  OATES  COVERED 


5.  FUMNNG  NUMBERS 
PE  602936 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AODRESS(ES) 

Naval  Undersea  Warfare  Center  Detachment 
39  Smith  Street 

New  London.  Connecticut  06320-5594 


8.  PERFORMING  ORGANIZAHON 
REPORTNUMBER 

TR  10,405 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Naval  Sea  Systems  Command 
2531  Jefferson  David  Highway 
Arlington,  VA  22242-5160 


10.  SPONSORING/MONITORING 
AGENCY  RS>ORT  NUMBER 


1 2a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximutn  200  words) 

In  the  deep  ocean,  long  range  active  sonars  are  often  volume  reverberation  limited  at  night  due  to  scattering 
from  the  gas-filled  swimbladders  of  mesopelagic  fish  in  near-surface  waters.  To  predict  the  detection  performance 
of  these  sonars  with  a  general  sonar  model,  volume  scattering  strength  profiles  ^z,f)  are  needed  as  inputs.  A 
general  volume  scattering  strength  model  (VSSM)  for  nighttime/daytime  Sv(z,f)  profiles  in  a  deep  water  faunal 
province  is  described.  Explicit  expressions  are  devetof^  tor  the  average  lerigth-depth  distrtoution  of  bladdered 
fish,  the  corresponding  swimbladder  radius-depth  distribution,  and  volume  scattering  strength  Syfz.f).  A  specific 
nighttime  near-surface  VSSM  is  developed  for  the  North  Sargasso  Sea  faunal  province  using  data  from  Ocean  Acre 
experiment  12.  A  VSSM  of  very  simple  form  is  found  capable  of  describing  measured  nighttime  profiles  at  3.85 
and  15.5  kHz  to  a  depth  of  225  m.  This  agreement,  however,  could  only  be  obtained  by  modifying  model 
parameters  obtained  from  trawl  data  to  correct  for  sampling  deficiencies,  principally  net  avoidance  by  the  larger  fish. 
Predictions  with  the  modified  VSSM  at  1  kHz  suggest  that  in  the  nighttime  0-50  m  surface  layer,  Sv*tovels  were  as 
high  as  -75  dB  due  to  scattering  from  these  larger  fish. 


14.  SUBJECT  TERMS 

Volume  scattering  strength  model  Ocean  Acre  scattering  layers  SwimUadder  target  strength 

Bioacoustic  volume  scattering  model  Mesopelagic  bladdered  fish  scatterers  Swimbladder  Q-faictor 


IS.  NUMBER  OF  PAGES 
46 


16.  PRICE  CODE 


17.  SECURITY  CLASSIRCATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  UkflTATON  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified  SAR 


NSN  7540-01 


PiwerOaS  by  ANSI  SU  Z3S-1S 
2M-102 


Table  of  Contents 


Page 


LIST  OF  ILLUSTRATIONS  ii 

I .  INTRODUCTION  1 

II.  A  PARAMETRIC  PHYSICAL  MODEL  FOR  PREDICTING  VOLUME  3 

SCATTERING  STRENGTH  IN  THE  DEEP  OCEAN  (VSSM) 

2.1  Biological  Modeling  Methodology  4 

2.2  Bioacoustic  Modeling  Methodology  10 

III.  A  SUMMER  NIGHTTIME  VSSM  FOR  NEAR-SURFACE  HATERS  OF  13 

THE  NORTH  SARGASSO  SEA  APPLICABLE  FOR  LOW  SONAR 

SEARCH  FREQUENCIES 

3.1  Biological  Conditions  at  Ocean  Acre  and  the  16 

Selection  of  a  VSSM 

3.2  Predictions  with  the  Single  Scattering  Layer  21 

Component  VSSM 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS  34 

V.  REFERENCES  36 

APPENDIX  A:  Theoretical  Representation  for  the  Al 

Characteristic  Pish  Length-Depth 
Density  Distribution  of  a  Scattering 
Layer  Component  in  the  VSSM 

APPENDIX  B:  Summary  of  the  Simplified  Andreeva  Bl 

Swimbladder  Target  Strength  Model 


1 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Representative  theoretical  Q(*^giean' 

Q(z)]^ax  curves  based  on  a  siaipli£ied  fora  of 
Andreeva's  (1964,  1972)  nodel  that  includes  only 
radiation  and  viscous  damping  losses.  Respective 
values  of  used  to  generate  these  curves 

are  shown  in  the  header.  Experiaental  data  from 
Table  1  are  included . 15 

2.  Resonance  frequency  vs  depth  domains  for  scatterers 

in  the  nominal  life  stage  size  classes  of  the 

Ocean  Acre  18  predominant  bladdered  species 

group . 19 

3.  TBSM  predictions  of  nighttime  S^(z,f (lower 

curve)  and  curve)  compared  with 


OA12N  experimental  profiles  of  Fisch  and  Dullea 
(1973)  taken  approximately  24  hours  apart:  Fig. 

(3a),  f-3.85  kHz;  Fig.  (3b),  f-15.5  kHz . 24 

4.  TBSM  predictions  of  nighttime 

Sv(z,f)Qmax  f-0.5  kHz  (Fig.  3a)  and  f-1.0  kHz 
(Fig.  3b) . 26 


5.  Number  density  distribution  n(z)  for  the  TBSM; 

Modified  TBSM.  The  characteristic  density 
distribution  n(z)  is  the  same  in  each  model  but 
the  respective  standing  crops  are  different: 

N-0.72/m^  and  N-4.93/m^ . 29 

6.  (a)  Fish  length  pdfs  p(L);  p(L)  for  the  TBSM; 

Modified  TBSM.  (b)  Corresponding  swimbladder 

radius  pdfs  g(R);  q(R)  for  the  two  models . 30 

7.  Modified  TBSM  predictions  f (lower  curve) 

and  ( 2 , f ) (upper  curve)  compared  with  OA12N 


ii 


Figure 


Page 


nighttine  experinental  profiles  of  Fisch  and 
Dullea  (1973)  taken  approximately  24  hours  apart: 

Fig.  (7a),  f«3.85  kHz;  Fig.  (7b),  f-15.5  kHz . 31 

8.  Predictions  with  the  Modified  TBSM  of  nighttime 
S^(z,f)Q^jjj  (lower  curve)  and 
curve)  for  f«0.5  kHz  (Fig.  8a)  and  f«1.0  kHz 
(Fig.  8b) . 33 


LIST  OF  TABLES 


Table  Page 

1.  Swimbladder  resonance  frequency  and  Q  measured  for 

various  live  swimbladdered  fish  with  swimbladders 
believed  to  be  intact . 14 

2.  Parameters  and  statistical  properties  of  the 

lognormal  pdfs  used  in  the  Trawl-Based  S^(z,f) 

Model  (TBSM)  and  the  Modified  TBSM . 23 


Accesion  For 

NTIS 

OTIC 

Unann 

Justific 

CRA&I 

TAB  O 

Dunced  □ 

ation 

•  w 

>'’i 

m  o 

ution  / 

Availability  Codes 

Dist 

Avail  c 
Spe 

r  d  /  or 
cial 

iii 


A  PREDICTIVE  PARAMETRIC  MODEL  FOR  VOLUME  SCATTERING  STRENGTH  IN 
THE  DEEP  OCEAN  WITH  AN  APPLICATION  TO  NIGHTTIME  NEAR-SURFACE 
SCATTERING  IN  THE  NORTH  SARGASSO  SEA  BELOW  16  kHz 

I .  INTRODUCTION 

In  the  deep  ocean,  volume  reverberation  below  20  kHz  is 
caused  principally  by  resonant  and  nonresonant  scattering  from 
the  gas-filled  swimbladders  of  numerous  species  of  mesopelagic 
fish  most  of  which  are  diel  vertical  migrators.  At  night,  a 
significant  fraction  of  these  bladdered  fish  scatterers  rise  from 
daytime  depths  to  feed  on  zooplankton  in  a  200  m  subsurface 
layer.  Under  commonly  occurring  propagation  conditions,  these 
near-surface  scattering  layers  can  significantly  lower  the 
detection  performance  of  long  range  2-5  kHz  search  sonars.  To 
predict  the  nighttime  detection  performance  of  such  sonars  for 
shallow  targets,  volume  scattering  strength  profiles  S^(z,f)  for 
these  layers  must  be  provided  as  inputs  to  sonar  performance 
prediction  models.  The  objective  is  to  show  how  S^-profiles  can 
be  parametrically  modeled  in  the  deep  ocean  and,  by  way  of 
illustration,  develop  a  nighttime  summer  model  for  the  North 
Sargasso  Sea. 

A  predictive  theoretical  model  for  S(z,f)  requires  submodels 
for  the  length-depth  distribution  of  the  scatterers,  the 
swimbladder  sizes  they  contain,  and  the  scattering  properties  of 
the  swimbladders.  First  we  present  a  general  theoretical  model 
for  predicting  S^(z,f)  in  a  deep  water  faunal  province  (Backus 
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and  Craddock  1977).  This  modal  is  than  spacifically  tailorad  to 
dascriba  nighttima  naar-surfaca  voluma  scattaring  in  tha  North 
Sargasso  Saa  during  summar.  Bioacoustic  data  from  OA12,  ona  of 
tha  most  succassful  of  tha  Ocaan  Acra  axparimants  (Brown  and 
Brooks,  1974),  ara  used  to  formulate  tha  modal.  Wa  note  hare 
that  the  Ocaan  Acra  measurement  site  is  a  1-dagraa  quadrangle 
southeast  of  Bermuda  centered  at  32*N,  64"W  and  that  it  is  a 
"typical”  location  in  the  North  Sargasso  Sea  faunal  province 
(Brooks  1972,  Backus  and  Craddock  1977). 

The  motivation  for  our  parametric  modeling  approach  stems  in 
part  from  the  limitations  of  discrete-depth  biological  sampling. 
In  the  Ocean  Acre  experiments  (Brown  and  Brooks  1974,  Gibbs  et 
al.  1987),  discrete-depth  biological  samples  were  obtained  with  a 
3-m  Isaacs  Kidd  Kidwater  Trawl  (IKNT)  and  concurrent  measurements 
of  S(z,f)  for  frequencies  between  3.85  and  15.5  kHz  were  obtained 
with  narrow-beam  transducers  (Fisch  and  Dullea  1973  and  Fisch 
1977).  While  the  seasonal  biological  properties  of  bladdered 
fish  at  the  Ocean  Acre  site  have  been  qualitatively  established 
by  these  trawl  measurements,  the  data  are  inadequate  for 
bioacoustic  prediction,  as  found  by  Brown  and  Brooks  (1974)  and 
Brooks  and  Brown  (1977).  This  is  due  to  the  sampling 
deficiencies  of  the  IKHT  caused  by  net  avoidance  and  escapement, 
further  complicated  by  patchiness  and  practical  limitations  on 
the  volume  of  water  that  can  be  sampled  at  each  depth. 

As  will  be  demonstrated,  a  "Trawl  Based  S(z,f)  Model"  (TBSN) 
with  parameters  obtained  directly  from  OA12N  data  leads  to 
unsatisfactory  agreement  between  predicted  and  measured  values  of 


2 


S^(z,f)  at  3.85  kHz,  as  night  be  expactad  for  tha  raasona  givan 
abova.  Bacausa  of  tha  paranatric  natura  of  tha  VSSM,  howavar, 
thoza  paranatara  dafining  tha  thaoratical  langth-dapth 
distribution  of  tha  scattarars  can  ba  adjustad  to  conpansata  for 
sampling  daficiancias.  Tha  "Modifiad  TBSN"  thus  obtainad  will  ba 
shown  to  produca  dramatically  battar  agraamant  batwaan  pradictad 
and  maasurad  S^-profilas  at  3.85  kHz  and  15.5  kHz.  S^-pradictions 
ara  also  mada  at  tha  nominal  fraquancias  0.5  and  1.0  kHz  at  which 
scattaring  from  masopalagic  fish  is  commonly  assumed  to  ba 
nagligibla.  The  large  differences  between  corresponding 
predictions  with  tha  two  models  at  3.85,  15.5,  0.5,  and  1.0  kHz 
ara  physically  explained  with  tha  help  of  a  resonance 
f reguency/depth  diagram  for  the  various  life  stage  size  classes 
of  the  bladdered  fish  population. 

II.  A  PARAHETRIC  PHYSICAL  MODEL  FOR  PREDICTING  VOLUME  SCATTERING 

STRENGTH  IN  THE  DEEP  OCEAN  (VSSM) 

In  section  2.1  we  first  formulate  an  idealized  theoretical 
model  for  the  length-depth  distribution  n(L,z)  of  swimbladdered 
fish.  The  approach  was  suggested  by  the  seasonal  depth 
distribution  patterns  of  bladdered  fish  observed  in  the  Ocean 
Acre  experiments  which  are  believed  to  be  fairly  typical  of  many 
deep  ocean  areas.  For  bioacoustic  modeling,  the  swimbladder 
radius-depth  distribution  n(R,z)  is  required,  and  the  procedure 
for  transforming  n(L,z)  to  n(R,z)  is  considered  next.  Then,  in 
section  2.2,  a  general  theoretical  formulation  for  the  volume 
backscatter  coefficient  s^(z,f)  is  given.  It  contains  a 
simplified  form  of  Andreeva's  (1964,  1972)  model  for  swimbladder 
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target  atrangth  in  which  contributions  of  hast  loss  to  daaping  at 
rasonanca  ara  ignorad.  Tha  haight  of  tha  swinbladdar  targat 
strangth  rasonanca  paak  is  datarainad  by  Qiz),  tha  invarsa 
damping  factor  at  rasonanca.  Raprasantativa  ainiaua,  aaan,  and 
aaxiaun  valuas  of  Q(s)  in  ganaral  accord  with  axpariaantal  valuas 
ara  proposad  for  S^pradictions  with  tha  VSSH.  Thasa  ara  obtainad 
from  tha  simplifiad  form  of  Andraeva's  formula  for  Q(z),  using 
suitabla  choicas  for  tha  complax  shear  modulus  of  fish  tissue. 

2.1  Biological  Modeling  Methodology 

Definitions  of  fish  lanqth-dapth  density  distributions 

and  associated  variables.  Consider  tha  swimbladdar  fish 

population  in  soma  deep  ocean  faunal  province  during  either 

nighttime  (N)  or  daytime  (D)  periods  when  the  distribution  of  the 

scatterers  is  approximately  steady  state  due  to  the  absence  of 

substantial  vertical  migration.  For  simplicityr  N/D  distribution 

conditions  are  assumed  to  be  spatially  and  temporally  homogeneous 

over  horizontal  x,y  planes  for  all  depths  of  interest.  We  define 

the  fish  length-depth  density  distribution  n(L,z)  such  that 

n(L,z)dLdz  represents  the  average  number  of  scatterers  with 

2 

lengths  L...L-<-dL  at  depths  z...z-«-dz  below  Im  of  sea  surface 
anywhere  within  the  faunal  province.  It  is  convenient  to  refer 

2 

this  distribution  to  the  average  number  of  scatterers  N  below  Im 
of  sea  surface,  the  "standing  crop,”  and  write 

(1)  n(L,z)  -  lik{h,z) 

where  n(L,z)  is  the  fish  length-depth  density  distribution  per 
unit  standing  crop.  We  assume  that  annual  variations  in  seasonal 
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distribution  proportios  sro  solely  due  to  annual  variations  in 
tha  seasonal  standing  crop  N.  This  pernits  us  to  interpret 
i^(L,s)  as  a  "characteristic”  seasonally  invariant  measure  of  the 
N/D  distribution  properties  of  the  faunal  province. 

The  quantity  {^(L,s)  may  also  be  interpreted  as  a  two 
dimensional  probability  density  distribution  function.  Prom 
probability  theory,  we  can  write  it  as 

(2)  A(L,s)  •  /i(z)p(L|z), 

where  ^(z)  is  a  probability  density  function  (pdf)  in  the  depth 
variable  z  and  p(L|z)  is  the  conditional  L-pdf  for  depth  z.  The 
marginal  L-pdf  associated  with  i^(L,z)  is 

m 

(3)  p(L)  -  J  A(L,z)dz  , 

o 

while  its  marginal  z-pdf  is 

«• 

(4)  (Hz)  -  J  k{h, z)6L  . 

o 

The  average  number  density  n(z),  i.e.,  the  average  number  of  fish 
per  unit  volume  at  depth  z,  is  given  by 

(5)  n(z)-iN^(z). 

Formulation  for  n(L,z)  in  the  VSSW.  Over  a  given  depth 
range  of  interest,  n(L,z)  is  expressed  as  a  sum  of  density 
distributions  nj^(L,z)  corresponding  to  i-1, 2, ..  .different 
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"scattering  layer  components" < 

(6)  n(LyS)  ■  n^CLdZ)  • 

A  scattering  layer  component  may  represent  an  individual 
bladdered  fish  species  or  a  group  of  such  species  treated  as  a 
single  "equivalent"  bladdered  species.  Distribution  properties 
of  mesopelagic  fish  are  size-dependent,  which  complicates 
biological  modeling.  Depending  on  the  frequencies  and  depth 
range  for  which  is  to  be  predicted,  not  all  size  classes  of 
fish  may  be  significant  contributors  to  scattering.  Restricting 
the  size  range  of  bladdered  fish  in  scattering  layer  components 
to  just  those  sizes  of  specific  acoustic  interest  may  simplify 
the  modeling  of  n^(L,z)  in  some  instances.  Seasonal  biological 
properties  of  the  bladdered  fish  fauna,  time-of-day,  and  the 
depths  and  frequency  ranges  for  which  S^(z,f)  is  to  be  predicted 
must  all  be  considered  in  chosing  the  appropriate  scattering 
layer  components  for  any  particular  VSSH. 

The  distribution  properties  of  the  i^^  scattering  layer 
component  are  formally  analogous  to  those  of  n(L,z)  described 
above.  Letting  be  the  standing  crop  and  nj^(L,z)«n^(z)p^(L|z) 
the  characteristic  distribution  of  the  i^^  scattering  layer 
component,  we  may  write  n(L,z)  as 

(7)  n(L,z)  -  5^  Nj^/lj(z)pj^(L|z)  . 

The  particular  characteristic  distributions  nj^(L,z)  employed 
in  the  prototype  VSSH  are  briefly  described  as  follows  (see 
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Appendix  A  for  details).  Consider  a  positive  valued  variate  e>0 

which  is  lognormally  distributed.  The  lognormal  pdf  of  9, 

2 

designated  here  by  X(9;ti,o  ),  is  given  by 

(8)  X(e;/#,o^)  -  [ee/5F]“^  EXP  [-(2o^)“^(lne-A/)^]  , 

2 

where  n  and  o  are  the  parameters  of  the  distribution.  As  shown 
in  Appendix  A,  the  characteristic  distribution 
i^j^(L,z)>n^(z)Pj^(L|  z)  for  an  i^^  scattering  layer  component  is 
defined  by  five  statistical  parameters,  ' ’’i '^xi ' '^xi '  **1  ^ 
terms  of  these  parameters, 

(9)  A^{z)  -  x[z;//^,ejj 

(10)  p^(L|z)  -  >^(l*?^xi|C'‘^xi|  J 


where,  in  the  first  distribution  parameter  of  Pj^(L|z),  Eq.  (10), 
C(  z  )«ln(  z/d) ,  with  d«l  m.  When  Pj^^O,  some  statistical  properties 
of  fish  length  such  as  the  mean,  e.g.,  vary  systematically  with  z 
throughout  the  i^**  layer,  due  to  the  change  in  form  of  p^(L|z) 
with  depth.  (Equations  for  the  z-dependence  of  the  mode,  median, 
and  mean  of  Pj^(L|z)  are  given  in  Appendix  A.)  If  the  random 
variates  L  and  z  are  statistically  independent  for  the  i^^ 
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scattering  layer  component,  then  and  we  have  a  "simple  i^^ 

scattering  layer  component".  For  such  components,  length 
distribution  properties  are  independent  of  depth,  as  may  be  seen 
by  putting  p^«0  in  Eq.  (10).  Then  the  conditional  pdf  Pj^(L|z) 
reduces  to  the  marginal  pdf  ),  Eq.  (A2),  in 

Appendix  A.  This  length  pdf  and  all  statistical  properties 
associated  with  it  remain  the  same  at  all  depths  throughout  the 
i^^  simple  scattering  layer  component. 

If  the  scattering  layer  components  of  n(L,z)  are  all  simple 
layers,  then  n(L,z)  in  Eq.  (7)  becomes 

(11)  n(L,z)  -  ^  Nj/l^(z)p.(L) 

w— 1  2  2 

•  X(  z  j /u ,  cf £  ^  ^ i  ^  *^xi  ^  * 

The  characteristic  distribution  nj^(L,z)  for  each  i  simple  layer 
is  specified  by  four  parameters,  »^xi '  “^xi  ^  scale  of 

its  contribution  to  n(L,z)  is  set  by  the  layer  standing  crop 
factor  . 

Transformation  of  n(L,z)  to  n(R,z).  Each  bladdered 
mesopelagic  fish  of  length  L  is  assumed  to  contain  a  fully 
inflated  swimbladder  with  an  equivalent  spherical  radius  R. 
Replacing  L  by  R  in  Eqs.  (7)  and  (11)  for  n(L,z),  we  have  the 
corresponding  formal  expressions 


(12) 

n(R,z)  - 

N^(*l^(z)qj.  (R|z) 

and 

(13) 

n(R,z)  - 

1 

Ni/ii(z)qi(R) 
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for  n(R,s).  To  obtain  explicit  expressions  for  n(R,s)  in  these 
equations,  transforsiations  Pj^(L|z} — ^q^(R|s)  and  p^iL) — ^q^iR) 
Bust  be  performed. 

Due  to  swiabladder  size  variability  in  aesopelagic  fish,  it 
is  desirable  for  purposes  of  transformation  to  use  a  more 
realistic  swimbladder  size  allometry  model  than  that  afforded  by 
an  allometric  law  regression  curve  R«aL^.  Bivariate  normal 
models  are  available  for  Ocean  Acre  bladdered  species  and 
bladdered  species  groups  (Saenger  1989).  If  such  models  are  used 
as  a  basis  for  transformation,  then  it  can  be  shown  (Saenger 

1988)  that  q.(R|z)  is  a  lognormal  pdf  with  distribution 

^  2 
parameters  depending  on  the  parameters  /^^i  |  C  ‘^xi  |  C 

pj^(L|z),  Eq.  (10),  and  other  quantities: 


^yi|C  “  *  ^^xi|C 

Vl?;  "  ^  xiK  ®y|x  ' 

where  a  and  b  refer  to  the  allometric  law  curve  R-aL^  associated 
with  the  bivariate  model,  and  s^jj^  is  a  constant  obtained  from 
the  model  that,  simply  put,  expresses  the  spread  of  swimbladder 
size  about  the  regression  curve.  For  the  special  case  pj^-0. 


(14) 

with 

(15) 

and 

(16) 
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where  we  are  dealing  with  a  simple  scattering  layer,  the 
transformed  distribution  q|(R)  obtained  from  Pj^(L),  Eq.  (A2),  is 
also  lognormal,  and  is  given  by 


(17)  q^CR)  -  » 

with 

(18)  -  Ina  + 
and 


(19) 


"yi 


■  b 


2_2 


xi  ®y|x 


Swimbladder  variability  extends  the  frequency  bandwidth  of  the 
scattering  so,  if  possible,  it  should  be  taken  into  account  to 
reduce  prediction  errors  in  S^.  (See  Saenger  1988  for  discussion 
and  examples). 


2.2  Bioacoustic  Modeling  Methodology 


The  volume  backscatter  coefficient  s^(z,f).  Randomly 
dispersed  or  loosely  aggregated  bladdered  fish  are  assumed  so 
that  scattering  contributions  from  individual  swimbladders  may  be 
summed  incoherently.  Following  conventions  introduced  by  Urick 
(1983),  the  volume  backscatter  coefficient  s^(z,f)  may  be 
expressed  as 

(20)  s^(z,f)  -  n(z)T(z,f )/d^,  m“^ 

2 

where  T(z,f)  is  the  average  target  strength  in  m  of  swimbladders 
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of  all  size  found  at  depth  z  for  an  ensonifylng  frequency  f. 

2  2 

(Some  prefer  to  OMit  d  *1b  from  this  definition,  in  which  case 
the  resulting  volusie  backscatter  coefficient  B^(z,f)  has 
dimensions  d"^.)  From  probability  theory,  T(z,f)  is  computed 
from 

m 

(21)  T(z,f)  -  Jq(R|z)TQ(R,z,f )dR,  m^ 

o 

2 

where  TQ(R,z,f)  is  the  target  strength  in  m  of  a  gas-filled 
swimbladder  of  radius  R  cm  at  depth  z  for  frequency  f.  Combining 
Eqs.  (20)  and  (21)  and  remembering  that  n(z)q(R|z)-n(R,z)*Ni^(R,z) , 
s^  may  be  written  as 

(22)  s^(z,f)  -  nJA(R,z) [TQ(R,z,f)/d^]dR 

o 

-  N^^(z,f) 


where  s^(z,f)  in  m”^  is  the  volume  backscatter  coefficient  per 
unit  standing  crop.  This  expression  points  up  the  central 
importance  of  the  characteristic  swimbladder  radius-depth 
distribution  n(R,z)  as  a  biological  input.  It  follows  from  our 
biological  assumptions  that  s^(z,f)  is  an  invariant 
characteristic  measure  of  seasonal  volume  scattering  properties 
in  the  faunal  province. 

The  volume  scattering  contribution  s^^(z,f)  from  each  i^*' 
scattering  layer  component  has  a  form  similar  to  that  for  s^(z,f) 
in  Eq.  ( 22 ) ,  viz . , 
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Due  to  the  assumption  of  incoherent  scattering,  s^j^-contributions 
from  each  i^^  component  can  be  added  to  give  s^(z,f).  with  L  and 
R  in  cm  and  z  in  m,  the  formal  theoretical  expressions  for 
S^(z,f)  in  the  VSSN  are  summarized  as  follows: 

(24)  S^(z,f)  -  10  log  [s^(z,f)d^],  dB  re  1/m^ 

(25)  s^(z,f)  -  jp  Nj^^i(z,f),  m"^ 


(26) 

^^i(z.t)  -  J 

c 

*qj(R|z)  [TQ(R,z,f  )/d^]dR,  m"^ 
> 

(27) 

/i^(z)  - 

2  -1 

X(zj//^,e£  ),  m 

(28) 

q^(R|z)  -  X| 

J'  ^  * 

For  a  simple  i^^  scattering  layer  component,  q^(R|z)  in  Eqs.  (26) 
and  (28)  is  replaced  by  qj^(R)  given  in  Eq.  (17).  At  present, 
full  computer  implementation  of  these  equations  exists  only  for 
simple  scattering  layer  components.  In  numerically  calculating 
1^vi(Z/£)»  Elq.  (26),  for  such  a  component,  digitally  sampled, 
truncated,  and  renormalized  pdfs  ^^^(z)  and  q£(R)  are  used,  and 
the  integral  is  evaluated  using  the  trapezoid  rule. 


The  simplified  Andreeva  model  for  svimbladder  target 
strength.  In  the  VSSM,  a  simplified  form  of  Andreeva's  (1964, 
1972)  model  for  swimbladder  target  strength  TQ(R,z,f)  is  used  in 
which  only  radiation  and  viscous  damping  losses  are  retained  in 
the  theoretical  expression  for  Q(z),  the  reciprocal  damping 


constant  at  resonance.  The  equations  of  this  nodel  as  used  in 
the  VSSK  are  sumarised  in  Appendix  B. 

In  modeling  volume  scattering  strength  in  the  deep  ocean,  it 

is  not  possible  to  measure  the  many  variables  that  determine  the 

Q  of  fish  swimbladders.  Therefore,  in  the  context  of  the 

simplified  Andreeva  model,  we  propose  the  adoption  of  provisional 

phenomenological  values  for  the  complex  shear  modulus 

of  fish  tissue  that  give  values  of  Q  in  reasonable  agreement  with 

experiment.  Some  available  experimental  data  are  shown  in  Table 

1.  Measurements  of  Q(z)  made  near  the  surface  are  of  the  order 

of  3  to  5.  Assuming  such  results  are  applicable  to  mesopelagic 

6  2 

fish,  and  taking  /t«^*l(10  )  dyn/cm  as  recommended  by  Batzler  and 
Pickwell  (1970),  Q(z)  was  computed  from  Eq.  B(10)  using  different 
values  of  By  trial-and-error ,  the  "Representative  Minimum, 

Mean,  and  Maximum  Q(z)  curves  shown  in  figure  1  were  obtained. 
These  curves  encompass  many  of  the  Q(z)  values  in  Table  1,  and 
the  trends  of  the  theoretical  Q(z)  curves  with  depth  are  at  least 
consistent  with  the  observations  of  McCartney  and  Stubbs  (1971). 
VSSM  predictions  displayed  below  utilize  one  or  other  of  these 
Q(z)  curves,  as  noted. 

III.  A  SUMMER  NIGHTTIME  VSSM  FOR  NEAR-SURFACE  WATERS  OF  THE 
NORTH  SARGASSO  SEA  APPLICABLE  FOR  LOW  SONAR  SEARCH 
FREQUENCIES 

The  VSSM  provides  a  general  formalism  for  modeling  the 
length-depth  distribution  n(L,z)  of  biological  scatterers  in  a 
faunal  province.  The  objective  here  is  to  select  the  simplest 
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Table  1.  Svlfflbladder  resonance  frequency  and  Q  measured  for  various  live  svimbladdered  fish  vith  svimbladders 
believed  to  be  intact. 
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(1)  The  example  given  is  claimed  to  be  "typical”  of  Q  values  between  3  and  5  obtained  from  aeasurem«its  on 
auichovy  and  goldfish. 

(2)  Condition  of  bladder  uncertain  -  fish  escaped. 

(3)  For  this  single  specimen,  f»  and  Q  are  adapted  values  at  each  depth. 


Q(z) 


Plot  of  Q(z)  —  S  imp lified  Andreeva  Model 

MU-R,  dynes/cm*  *2  t.0000-f06  1  .000Cm>6  1.0000+06 

MU- 1 M.  dyne  s/am*  2  2. 5000+06  1.5000+06  1.0000+06 


Figure  1.  Representative  theoretical  Q(*)niin'  ®^*^mean'  ®^*^iiiax 

curves  based  on  a  simplified  form  of  Andreeva's  (1964^  1972) 
model  that  includes  only  radiation  and  viscous  damping 
losses.  Respective  values  of  used  to  generate 

these  curves  are  shown  in  the  header.  Experimental  data  from 
Table  1  are  included. 
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model  capable  of  predicting  S(z,£)  between  the  surface  and  a 
depth  of  the  order  of  200  m  for  sonar  frequencies  up  to  about  15 
kHz.  This  requires  consideration  of  the  broad  features  of  OA12N 
biological  data  and  their  physical  implications.  Unfortunately, 
average  trends  in  the  length-depth  distributions  of  bladdered 
fish  are  obscured  by  the  unknown  efficiency  of  the  IKMT  net  and 
sampling  fluctuations  accentuated  by  the  low  densities  of 
bladdered  fish  and  patchiness.  Consequently,  the  choice  of  model 
is  somewhat  subjective.  The  ultimate  justification  for  any 
particular  model  is  that  it  involve  as  few  theoretical  parameters 
as  possible  and  that  it  is  capable  of  predicting  experimental 
values  of  S(z,f).  In  the  section  below,  a  very  simple  model  for 
n(L,z)  is  proposed  based  on  distribution  properties  inferred  from 
OA12N  trawl  data.  Predictions  with  this  model  are  made  in  the 
following  section  where  it  is  shown  that,  due  primarily  to  net 
avoidance  by  the  larger  fish,  model  parameters  derived  from  trawl 
data  must  be  modified  to  obtain  realistic  predictions. 


3.1  Biological  Conditions  at  Ocean  Acre  and  the  Selection  of  a 

VSSM 

Using  trawl  data  from  the  major  Ocean  Acre  experiments 
together  with  swimbladder  size  allometry  information.  Brown  and 
Brooks  (1974)  identified  55  predominant  swimbladdered  fish  taxa 
at  Ocean  Acre.  The  size-depth  distribution  properties  of  these 
fish  during  the  spring  (OA10,14),  summer  (OA12),  and  winter 
(OA13)  were  studied  in  detail  by  Brooks  and  Saenger  (1991).  In 
the  OA12N  nighttime  0-200  m  layer,  these  fish  have  individual 
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overlapping  depth  niches  and  a  wide  range  of  fish  size  is  found 
at  each  depth.  This  configuration  is  believed  to  represent  a 
basic  ecological  strategy  for  reducing  competition  for 
zooplankton  food  resources  during  nighttime  feeding  periods. 
Because  of  the  small  numbers  caught,  length-depth  distributions 
of  individual  bladdered  fish  species  are  ill-defined 
experimentally  so  a  "species  group"  approach  to  bioacoustic 
modeling  is  necessary. 

Brooks  and  Saenger  (1991)  found  that  18  particular 
predominant  bladdered  species  taken  as  a  group  had  diel 
length-depth  distribution  properties  that  closely  represent  the 
distribution  properties  of  the  entire  night/day  catches  of 
swimbladdered  fish  at  Ocean  Acre  in  spring,  summer  and  winter. 
This  18  predominant  bladdered  species  group  (18PBSG)  was 
partitioned  into  larval-postlarval  (L/PL),  juvenile  (J),  subadult 
(SA),  and  adult  life-stage  size  classes.  In  0A12,  the  nighttime 
composition  of  the  18PBSG  group  found  between  the  surface  and  a 
depth  of  225  m  was  L/PL(40.5%),  J(43.5%),  SA(14.5%),  and  A(1.5%). 
The  characteristic  density  distribution  n  ( z )  for  each  life  stage 
size  class,  normalized  with  respect  to  the  0-225  m  depth 
interval,  is  unimodal  and  the  distributions  for  all  size  classes 
are  fairly  similar.  Also,  the  experimental  fish  length 
distributions  p(L|Z|^)  at  the  sample  depths  Z|^«25,  125,  175,  and 
225  m  are  unimodal  and  do  not  appear  to  depend  strongly  on  depth. 
(No  18PBSG  fish  were  taken  at  Zj^-O.)  On  the  basis  of  this 
subjective  identification  of  major  trends  in  the  trawl  data,  it 
is  proposed  that  n(L,z)  for  the  18PBSG  be  modeled  as  a  simple 
scattering  layer  taking 
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(29) 


n(L,s)  -  nA(x)p(L) 


Further,  to  represent  the  swiabladder  else  alloaetry  properties 

of  this  group,  it  seeas  most  appropriate  to  use  an  "equivalent 

bivariate  normal”  model  developed  by  Saenger  (1969)  for  a 

subgroup  of  50  Ocean  Acre  nonregressive  swimbladdered  fish  taxa 

drawn  from  the  list  of  55  OA  predominant  swimbladdered  fish  taxa 

given  by  Brown  and  Brooks  (1974).  Regression  equation  constants 

(a,b)  and  s^.  from  this  model  are  found  at  the  bottom  of  Table  2 

below.  The  swimbladder  radius-depth  distribution  corresponding 

to  n(L,z)  in  Eq.  (29)  is  n(R,z)<-NA(z)q(R) ,  and  the  parameters 
2 

(//y,ey  )  of  the  lognormal  pdf  q(R)  are  obtained  from  those  of 
p(L)  using  Eqs.  (18)  and  (19).  Finally,  the  expression  for 
s^(z,f)  in  this  VSSH  is  obtained  by  substituting  Nfi(z)q(R)  for 
Nfi(R,z)  in  Eq.  (22). 

Physical  implications  and  possible  limitations  of  this  very 
simple  model  are  conveniently  discussed  with  the  help  of  figure 
2,  a  resonance  frequency/depth  diagram  for  the  different  size 
classes.  The  9  cm  representative  upper  bound  for  the  A  class  is 
suggested  by  modeling  results  in  the  next  section.  Nominal 
L-boundary  values  between  each  size  class  were  transformed  to 
corresponding  R-values  using  R-aL^,  and  resonance  frequencies 
fQ(R,z)  were  then  plotted  as  a  function  of  depth  using  Andreeva's 
formula,  Eq.  (B5).  From  Eq.  (29)  it  is  seen  that  Ni^(z) 
represents  the  average  number  of  bladdered  fish  of  all  sizes  at 
depth  z.  The  relative  proportions  of  the  different  size  classes 
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res ,  kHz 


Swimblodder  Resonance  Frequency  vs  z  for  Various  Fish  Lengths  L 
L ( cm)  1 . 3500+00  2 . 2500+00  3  3500+00  9 . 0000+00 

Porometers:  o  =  3.3840-02  b  =  8.9190-01  MU— R  =  1.0000+06  dynes/cm* *2 


Z  .  m 


Figure  2.  Resonance  frequency  vs  depth  domains  for  scatterers  in  the 
nominal  life  stage  size  classes  of  the  Ocean  Acre  18 
predominant  bladdered  species  group. 
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are  detemined  by  p(L)  and  are  the  sane  at  each  depth  in  this 
model.  For  a  given  sonar  frequency  f,  contributions  to  S^(z,£) 
from  fish  of  a  particular  size  will  be  greatest  at  those  depths 
where  the  resonance  frequency  f^CR^z)  of  their  swimbladders  is 
equal  to,  or  near,  f.  Since  TQ(R,z,fQ)  «  from  Eq.  (B3), 

such  contributions  from  large  fish  at  low  kHz  frequencies  can  be 
very  important,  even  when  these  fish  are  present  in  relatively 
low  numbers.  Therefore  the  accuracy  of  low  frequency  predictions 
with  this  VSSN  depends  critically  on  the  tail  portion  of  the 
model  pdf  p(L). 

The  size  pass  band  limitations  of  the  IKNT  and  the  expected 
large  summer  increase  of  L/PL  fish  at  Ocean  Acre  must  both  be 
considered  in  assessing  the  probable  effectiveness  of  the 
proposed  model  for  n(L,z).  Gibbs  and  Karnella  (1987)  tentatively 
estimated  the  range  of  fish  size  efficiently  sampled  by  the  3-m 
IKHT  to  be  1-5  cm.  With  regard  to  the  larger  fish,  comparison  of 
0A12N  non-discrete  depth  oblique  catches  with  the  much  larger 
Engel  trawl  suggest  that  adult  fish  are  inadequately  represented 
in  the  IKHT  catches.  Brooks  and  Brown  (1978)  reported  that  while 
the  range  of  fish  size  obtained  with  both  nets  was  similar,  the 
mean  lengths  of  fish  caught  with  the  Engel  trawl  were  commonly 
twice  as  large.  In  particular  (private  communication),  they 
found  the  nighttime  size  distribution  for  a  group  of  6  species 
(which  incidentally  happen  to  belong  to  the  18PBSG)  to  be 
unimodal,  with  the  mode  of  the  IKHT  distribution  at  1.2  cm  and 
that  for  the  Engel  distribution  at  2.6  cm.  Thus,  it  is  very 
probable  that  the  larger  adult  fish  are  not  present  in  the  IKHT 
catch.  Also,  it  is  very  probable  that  the  L/PL  size  class  is 
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undersanpled,  since  the  estimated  decrease  in  catch  efficiency  of 
the  IKMT  occurs  over  much  of  the  L/PL  size  range. 

If  the  summertime  L/PL  size  class  is  actually  a  larger 
fraction  of  the  18PBSG  than  is  indicated  by  the  IKHT  catch,  it 
might  be  necessary  to  use  a  more  general  n(L,z)  model  than  that 
proposed  in  Eg.  29.  In  view  of  the  size-dependent  differences  in 
life  style  expected  of  the  L/PL  and  (J-t-SA-fA)  size  groups,  one 
might  treat  these  size  groups  as  independent  scattering  layer 
components.  Or,  individual  size  classes  might  be  treated  as 
independent  scattering  layer  components  in  an  even  more  refined 
model.  However,  trawl  data  provide  no  definitive  evidence  for 
the  need  of  one  or  other  of  these  more  complex  models,  so  we 
shall  explore  the  predictive  capability  of  the  single  scattering 
layer  component  model,  the  simplest  of  the  various  possible 
VSSHs. 

3.2  Predictions  with  the  Single  Scattering  Layer  Component  VSSN 

The  Trawl-Based  S  Model  (TBSM).  Predictions  in  the  TBSM  are 

made  with  model  parameters  derived  solely  from  0A12N  IKMT  sample 

data  for  the  18PBSG  taken  at  the  depths  25,  125,  175,  and 

225  m.  No  attempt  has  been  made  to  correct  the  parameters  for 

2  2 

sampling  deficiencies.  [N,/e/,o  ^  calculated  as 

_  -r 

follows.  Estimates  of  N  and  the  moments  z  and  z  of  the 

experimental  normalized  density  distribution  referred  to  the 

0-225  m  depth  interval  were  obtained  using  procedures  described 

2  A 

in  Brooks  and  Saenger  1991.  The  parameters  (/u,a  )  of  n(z)  were 

_  ~T 

then  derived  from  z  and  z  by  the  method  of  moments  (Aitchison 
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2 

and  Brown  1957).  The  parameters  )  of  p(L)  were  Tsompu ted 

2 

as  density-weighted  averages  of  latter 

were  calculated  by  applying  the  method  of  moments  to  experimental 

—  “T  2 

values  of  L(Z|()  and  L  (Sj^).  Then  the  parameters  (f/ytOy  ) 

defining  q(R)  were  obtained  from  using  Eqs.  (18)  and 

(19)  with  allometric  constants  a,b,  and  from  the  50 

swimbladdered  fish  taxa  bivariate  normal  model  mentioned  above. 

Table  2  contains  the  TBSH  parameters  so  obtained  and  some 

statistical  properties  of  the  distributions  (^(z)  and  p(L). 

TBSH  predictions  S(z,f and  S ( z , f ) for  3.85  kHz  are 
compared  with  OA12N  experimental  S^-profiles  obtained  on  two 
successive  nights  by  Fisch  and  Dullea  (1973)  in  Fig.  (3a).  The 
agreement  is  very  poor.  From  our  previous  considerations,  some 
of  the  shortfall  in  overall  level  may  be  supposed  due  to  the 
model  value  for  standing  crop  N  being  too  low.  However,  the 
increasing  divergence  of  the  depth  trend  of  predicted  and 
experimental  values  is  consistent  with  the  expectation  that 
important  scattering  contributions  from  larger  fish  are  missing 
from  the  model.  As  seen  from  figure  2,  contributions  from  adult 
fish  are  expected  to  be  sizable  at  3.85  kHz,  not  only  because  of 
the  large  size  of  their  swimbladders,  but  because  scattering 
occurs  at  a  frequency  equal  to,  or  not  too  far  from,  the 
resonance  frequency  of  their  swimbladders. 

At  15.5  kHz,  Fig.  (3b),  agreement  with  experiment  is  still 
poor,  but  there  is  now  less  discrepancy  in  overall  level  and 
depth  trends  of  predicted  and  experimental  values.  From  figure 
2,  this  may  be  tentatively  explained  by  noting  that  at  this 
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Table  2.  Paraaeters  and  statistical  properties  of  the  lognoraal  pdfs  used  in  the  Travl>Based  S^(s»f)  Model 
(TBSM)  and  the  Modified  TBSH.* 
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(see  text). 

Eqs.  (18)  and  (19)  vere  used  to  transfom  p(L)  to  q(R)  and  p(L)  to  q(R)  vith  a-3.384(10~ 
b-0.8919,  and  s^,  >0.1207  (Ocean  Acre  50  nonregressive  bladdered  fish  taxa  uodel). 


0A12N  TRAWL-BASED  S(z.f)  PREDICTION  MODEL.  f=3.85  kHz 
COMPARISON  OF  EXPERIMENTAL  WITH  MINIMUM  AND  MAXIMUM  PREDICTED  VALUES 


OA12N  TRAWL -BASED  S(z,f)  PREDICTION  MODEL.  f=15.5  kHz 
COMPARISON  OF  EXPERIMENTAL  WITH  MINIMUM  AND  MAXIMUM  PREDICTED  VALUES 


Depth  z .  m 


Figure  3.  TBSM  predictions  of  nighttime  z ,  f  (lower  curve)  and 

(  z , f ) Qjugjjj  (upper  curve)  compared  with  OA12N  experimental 
profiles  of  Fisch  and  Dullea  (1973)  taken  approximately  24 
hours  apart:  Fig.  (3a),  f«3.85  kHz;  Fig.  (3b),  f“15.5  kHz. 
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higher  frequency,  and  for  the  depths  where  comparison  is 
possible,  dominant  scattering  contributions  very  likely  came  from 
J  and  SA  fish  that  on  the  whole  are  expected  to  have  been  sampled 
more  efficiently. 

TBSM  predictions  for  f>0.5  and  1.0  kHz  are  sho%im  in  figures 

(4a)  and  (4b),  respectively.  Xn  each  plot,  the  two  curves  are 

practically  identical,  since  Q  does  not  affect  the  form  of 

S^(z,f)  at  these  low  frequencies.  Here  Rayleigh-type  scattering 

2  -4  4  2 

is  occurring  so,  in  m  ,  TQ(R,z,f )s>10  (f/f^)  R  ,  with  R  in  cm. 

This  explains  the  systematic  12  db  difference  in  predicted  values 
of  S^(z,f)  at  the  two  frequencies.  The  lowness  of  the  predicted 
levels  is  presumably  caused  by  the  absence  of  scattering 
contributions  from  the  larger  fish  that  would  be  the  dominant 
contributors  to  at  these  low  frequencies. 

A  Modified  TBSM  with  model  parameters  corrected  for  net 
sampling  deficiencies.  The  choice  of  TBSH  parameters  requiring 
modification  to  obtain  agreement  of  VSSN  predictions  with 
experiment  involves  a  number  of  physical  considerations.  First, 
it  seems  reasonable  to  assume  that  n(z)  obtained  from  the  TBSH 
provides  a  reasonable  approximation  for  the  actual  characteristic 
density  distribution.  This  assumption  is  prompted  by  the 
following  observations:  1)  ^(z)  can  be  expressed  as  a  sum  of 
contributions  (N-/N)ft  (z)  from  each  s^^  size  class,  where  N  is 

So  S 

the  standing  crop  of  the  s^**  size  class;  2)  n  (z)-n  (z)/N  ,  by 

its  definition,  tends  to  be  independent  of  a  net  sampling 

efficiency  factor  for  the  s^**  size  class;  3)  ft_(z)  for  the  L/PL 

s 

and  J  size  classes  are  virtually  identical;  and  4)  these  two  size 
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RANGE  OF  OA12N  TRAWL -BASED  S(z.f)  PREDICTIONS  FOR  f»0.5  kHz 


Figure 


RANGE  OF  OA12N  TRAWL-BASED  S(2.f)  PREDICTIONS  FOR  f=1 


.  TBSM  predictions  of  nighttine  and  z , f 

f-0.5  kHz  (Fig.  3a)  and  f-1.0  kHz  (Fig.  3b). 
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classes  constitute  84%  of  the  18PBSG  catch  between  the  surface 

2 

and  225  m  depth.  Assuming,  then,  that  in, a  )  are  adequate,  only 

2 

the  parameters  '  need  modification,  in  estimating  these 

parameters,  we  primarily  utilize  experimental  S^-data  at  3.85 
kHz,  the  lowest  frequency  for  which  experimental  data  from  0A12N 
are  available,  since  scattering  contributions  from  larger  fish 
missed  by  the  IKMT  are  best  represented  in  this  data. 

The  average  length-depth  density  in  the  Modified  TBSM  is 
designated  by 


(30)  n(L,z)  -  M(z)p(L) 

The  parameters  of  p(L)  were  obtained  by  systematically  modifying 
those  of  p(L)  in  such  a  way  as  to  increase  the  relative 
proportion  of  larger  fish  in  the  distribution  without  unduly 
suppressing  the  smaller  size  class  components.  This  was  done  by 
successively  increasing  the  3  cm  97.5%  fractile  point  of  p(L)  to 
6  cm,  9  cm,  15  cm...  The  corresponding  transformed  pdfs  q(R) 
were  used  to  predict  3*85  and  15.5  kHz.  For  each 

case,  the  standing  crop  parameter  N  was  obtained  by  equating 
interval  column  strengths  for  the  23-221  m  depth  interval 
computed  from  the  predicted  and  experimental  3.85  kHz 
S^-profiles.  Selection  of  the  "optimum"  Modified  TBSM  was 
necessarily  somewhat  subjective.  Properties  satisfied  by  the 
model  finally  chosen  were  the  following: 
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1.  Physical  reasonableness  of  the  97.5%  fractile  points  of 

•  • 

the  p(L)  and  q(R)  distributions; 

2.  Acceptable  agreement  in  the  depth  trends  of  3.85  kHz 
^v^*'^^Qmean  experimental  S^-profiles,  apart  from 
absolute  levels; 

3.  N  such  as  to  give  better  agreement  in  predicted  levels 
at  15.5  kHz  than  were  obtained  with  the  TBSH. 

Model  constants  and  statistical  information  for  the  Modified  TBSM 
are  given  in  Table  2. 

The  number  density  distributions  n(z)  in  the  TBSM  and 
Modified  TBSM  are  shown  in  figure  5.  It  is  seen  that,  to  get  the 
agreement  with  measured  S^-values  shown  below,  a  larger  standing 
crop  is  required.  Also,  the  relative  proportion  of  larger  fish 
must  be  increased,  as  seen  in  figure  (6a).  These  results  imply 
that  both  the  absolute  values  of  fish  density  and  the  size  ranges 
of  fish  caught  with  the  IKMT  are  seriously  biased.  Corresponding 
swimbladder  size  pdfs  are  compared  in  figure  (6b). 

The  relative  increase  in  large  swimbladder  radii  in  the 
Modified  TBSM  has  a  substantial  impact  on  the  predicted  levels  of 
at  3.85  and  15.5  kHz,  as  may  be  seen  by  comparing  Figs.  (3a) 
and  (7a)  and  Figs.  (3b)  and  (7b).  These  results  indicate  that 
the  single  scattering  component  VSSM  can  provide  acceptable  fits 
to  measured  S^-profiles  at  frequencies  where  different  size 
classes  in  the  bladdered  fish  population  are  the  dominant 
scatterers.  This  implies  that  the  length-depth  distribution  of 
bladdered  fish  given  by  the  Modified  TBSM  is  very  likely  fairly 
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0A12N  NUMBER  density  FOR  THE  TBSM  AND  MODIFIED  TBSM 
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Figure  5.  Number  density  distribution  n(z)  for  the  TBSM;  Modified  TBSM. 

The  characteristic  density  distribution  n(z)  is  the  same  in 
each  model  but  the  respective  standing  crops  are  different: 
N-0.72/m^  and  N»4.93/m  . 
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Figure  6 


0.1  0.2  0.3  0.4 

Swimbladder  Radius  R,  cm 

(a)  Fish  length  pdfs  p(L);  p(L)  for  the  TBSM;  Modified  TBSM. 

(b)  Corresponding  swimbladder  radius  pdfs  q(R);  q(R)  for  the 
two  models. 


0A12N  modified  trawl -based  S{z.1)  PREDICTION  MODEL,  f-3.85  kHz 
COMPARISON  OF  EXPERIMENTAL  WITH  MINIMUM  AND  MAXIMUM  PREDICTED  VALUES 
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0A12N  modified  trawl -based  S(z.f)  PREDICTION  MODEL,f=15.5  kHz 
COMPARISON  OF  EXPERIMENTAL  WITH  MINIMUM  AND  MAXIMUM  PREDICTED  VALUES 
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Figure  7.  Modified  TBSM  predictions  S^(z,f (lower  curve)  and 
S^(z,f (upper  curve)  compared  with  OA12N  nighttime 
experimental  profiles  of  Fisch  and  Dullea  (1973)  taken 
approximately  24  hours  apart:  Fig.  (7a),  f«3.85  kHz; 
Fig.  (7b),  f-15.5  kHz. 


realistic,  and  that  substantially  siore  fish,  particularly  larger 
ones,  are  present  than  is  indicated  by  discrete  depth 
measurements  with  the  IKNT. 

Finally,  predictions  of  S^(s,f)  at  0.5  and  1.0  kHs  are  shown 
in  figures  (8a)  and  (8b).  Comparison  with  the  corresponding  IBSM 
predictions  in  figures  (4a)  and  (4b)  shows  a  marked  increase  in 
predicted  levels  at  both  frequencies  and  also  changes  in  profile 
shape,  particularly  in  the  1.0  kHz  profile,  due  to  the  inclusion 
of  scattering  contributions  from  larger  fish.  Since  the 
resonance  frequencies  f^  of  the  swimbladders  of  these  larger  fish 
are  not  very  far  above  0.5  and  1.0  kHz  at  shallow  depths, 
differences  in  Q  now  have  a  perceptible  effect  on  predicted 
values  at  these  depths,  as  seen  in  figures  (8a)  and  (8b). 
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IV.  Conclusions  &  Recosuiendations 


In  deep  water,  long  range  active  sonars  can  be  volume 
reverberation  limited  at  night  due  to  scattering  from  small 
mesopelagic  swimbladdered  fish  present  in  near-surface  waters. 

The  effect  of  these  scatterers  on  detection  performance  in  a 
given  propagation  situation  can  be  predicted  with  a  general  sonar 
model  if  a  scattering  strength  profile  S^(z,f)  is  available. 

A  parametric  volume  scattering  strength  model  (VSSM)  for 
predicting  night time /day time  volume  scattering  strength  profiles 
S^(z,f)  in  a  deep  water  faunal  province  for  frequencies  below  20 
kHz  has  been  presented.  The  volume  backscatter  coefficient 
s^(z,f)  is  formulated  as  the  product  of  a  normalized  volume 
backscatter  coefficient  representing  the  characteristic 
scattering  properties  of  the  province  for  a  particular  season, 
and  the  seasonal  standing  crop  of  bladdered  fish  scatterers  which 
generally  varies  from  year-to-year.  The  similarity  of 
nighttime/daytime  characteristic  number  density  distributions  of 
bladdered  species  groups  in  different  ocean  locations  shown  by 
Brooks  and  Saenger  (1991)  and  the  general  similarity  of 
scattering  layers  in  the  deep  ocean  suggests  that  the  VSSH  may 
prove  useful  for  modeling  S^(z,f)  in  many  deep  water  areas. 

To  test  the  VSSH,  a  summer  nighttime  model  for  near-surface 
scattering  has  been  developed  for  the  North  Sargasso  Sea  faunal 
province  using  Ocean  Acre  12  IKNT  trawl  data  and  experimental 
S^(z,f)  profiles.  Results  obtained  are  notable  for  several 
reasons.  First,  the  simplest  form  of  VSSH  was  found  capable  of 
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describing  the  0-225  ■  portions  of  experimental  S^-pro£iles  at 
3.85  and  15.5  kHz.  This  could  only  be  done,  however,  by 
modifying  the  model  parameters  calculated  from  trawl  data  that 
describe  the  length-depth  distribution  of  scatterers  caught  with 
the  IKMT.  The  need  to  substantially  Increase  both  the  standing 
crop  and  the  relative  proportion  of  larger  individuals  in  the 
model  bladdered  fish  population  in  order  to  predict  measured 
values  of  at  3.85  and  15.5  kHz  further  confirms  the  long 
suspected  limitations  of  discrete  depth  sampling  with  the  3-m 
IKHT.  Second,  the  fact  that  dominant  scattering  contributions  at 
3.85  and  15.5  come  from  different  size  classes  of  the  population 
suggests  that  the  actual  OA12N  length-depth  distribution  of 
bladdered  fish  between  the  surface  and  225  m  is  adequately 
modeled  with  this  simple  model.  Finally,  model  predictions  for 
lower  frequencies  suggest  that,  due  to  relatively  strong 
scattering  contributions  from  the  larger  mesopelagic  fish, 
S^-levels  in  the  nighttime  0-50  m  layer  were  as  high  as  -75  dB  at 
1  kHz . 

It  is  recommended  that 

1.  The  generality  of  the  single  scattering  layer  component 
VSSM  for  nighttime  near-surface  scattering  be  tested  with 
bioacoustic  data  from  other  Ocean  Acre  experiments  and  available 
bioacoustic  data  from  other  deep  water  areas;  and 

2.  Computer  implementation  of  the  VSSN  be  completed  and 
documented  to  facilitate  generation  of  deep  water  model  S^(z,f) 
profiles  for  use  in  Navy  sonar  performance  prediction  models. 
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Appendix  A:  Theoretical  Representation  for  the  Characteristic 
Fish  Length>Depth  Density  Distribution  of  a 
Scattering  Layer  Conponent  in  the  vssn 

As  shown  in  Eq.  (7),  the  fish  length-depth  density 
distribution  n(L,z)  consists  of  a  sum  of  N^-weighted  characteristic 
density  distributions 

(Al)  A^(L,z)  -  A^(z)p^(L|z) 

for  scattering  layer  components  i-1,2,....  He  require  an 
explicit  parametric  model  for  n^(L,z).  In  the  present  form  of 
the  VSSN,  the  model  chosen  is  based  on  general  mathematical  and 
physical  considerations,  rather  than  specific  assumptions  as  to 
the  physical/statistical  response  of  the  scatterers  to  a 
particular  environmental  variable  or  variables.  Let  x«ln(L/l) 
and  C*ln(z/d)  be  logarithmic  measures  of  fish  length  and  depth, 
respectively,  with  1-1  cm  and  d-1  m.  Dropping  the  subscript  i  in 
what  follows  for  notational  convenience,  we  propose  the  following 

Fundamental  Hypothesis:  in  terms  of  the  length  variate  x 
and  the  depth  variate  C,  the  characteristic  density  distribution 
g(x,C)  of  a  scattering  layer  component  is  bivariate  normal  with 
parameters 

The  properties  of  the  bivariate  normal  density  distribution 
are  well  known.  (A  simple  exposition  sufficient  for  present 
purposes  may  be  found,  e.g.,  in  Cramer  1955.)  The  marginal  x  and 
C  pdfs  and  the  conditional  x  pdf  of  g(x,C)  are  normal,  so  the 
corresponding  pdfs  in  the  positive  variates  L  and  z  are  lognormal 


Al 


(Aitchison  and  Brown  1957): 


(A2) 

p(L) 

(A3) 

k{z) 

-  X^z;/i,e^j,  and 

(A4) 

P(L|Z) 

■  '‘('■’''xic-'xic)' 

(AS) 

^xK 

(A6) 

®x|i; 

For  a  general  scattering  layer  component  (p>*0),  the 
conditional  pdf  p(L|z)  and  its  statistical  properties  change 
systematically  with  depth  z«de^.  For  C-//,  z-z  -de^,  the  median 

A  “ 

of  n(z).  One  can  show  that 

<*®>  ''median'*’  -  (*/*«)"’*''” 

(A9)  ‘'mean'*’  ’  '  J]  (a/zJ''V». 

The  coefficient  of  variation  of  p(L|z)  (i.e.,  the  standard 

deviation  divided  by  the  mean)  has  the  constant  value 


(AlO) 


For  a  simple  scattering  layer  component  (p«0),  the 


A2 


conditional  pdf  p(L|z),  Eq.  (A4)/  reduces  to  the  narginal  pdf 
p(L),  Eq.  (A2),  and  all  statistical  properties  of  L  are 
independent  of  depth.  Formulas  for  the  mode/  median,  mean,  and 
coefficient  of  variation  of  p(L)  are  obtained  by  evaluating  the 
right  hand  sides  of  Eqs.  (A7)  -  (AlO),  respectively,  taking  p-0. 
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Appendix  B:  Suaaary  of  the  Simplified  Andreeva  Swimbladder  Target 
Strength  Model 

The  basic  equations  of  the  model  (Andreeva  1964,  1972)  are 
conveniently  summarised  as  follows:  Let 

(Bl)  G(z)  -  [l/(4R2)j [(3Yp(*)+4^^)/p^j  , 

where  Y*Cp/c^  for  the  gas  in  the  swimbladder,  P(z)  is  hydrostatic 
pressure  at  depth  z,  is  the  real  part  of  the  complex  shear 
modulus  of  fish  tissue,  and  is  the  density  of  fish 

tissue,  approximately  that  of  seawater.  Also,  letting 

(B2)  b(z)  -  1-[2q2(z)]"^  , 

we  may  write  the  target  strength  TQ(R,z,f)  of  a  gas-filled 
swimbladder  of  radius  R  at  depth  z  for  frequency  f  as 

(B3)  TQ(R,z,f)  -  RV[R^-2b(z)R5R^+Rj]  , 

where 

(B4)  R*(z,f)  -  /G(z)/f  , 

the  "resonance  radius".  (At  depth  z,  for  an  ensonifying 
frequency  f,  a  swimbladder  with  radius  R^  will  be  at  resonance.) 
We  note  here  that  the  "inverse"  of  this  equation  is 

(B5)  fo(R»z)  -  /g(z)/r  , 

Andreeva's  formula  for  the  resonance  frequency  fQ  of  a 
swimbladder  of  radius  R  at  depth  z. 

Neglecting  thermal  damping,  Andreeva's  model  for  Q(z)  may  be 


Bl 


written  as 


[0(»)]'‘-  ri,o  +  "vo 

where  the  danping  constant  at  resonance  due  to  radiation  losses 
is 

(B7)  hjo  •  "  <2sfQ/c)R 

-  (2r/c)  /G(s) ■  , 

using  (b5),  and  that  due  to  viscous  losses  is 

(B8)  -  4pj,[3tP(*)+4„J-1 

-  - 

using  (Bl). 


For  practical  calculations  with  the  VSSM,  we  use  the 

hydrostatic  formula  for  P(z)  and  the  values  P(0)«1.013( 10^)  dyn/cm 

3  2  3 

p^-1.025  g/cm  ,  g-980.6  cm/s  ,  y*1.40  for  air,  p^«1.08g/cm  , 

5  2 

c-1.521(10  )  cm/sec.  R  and  R^  are  in  cm,  f  in  kHz,  Tq  in  cm  , 

2 

and  are  in  dyn/cm  .  (In  the  VSSN,  Tq  from  Eq.  (B3)  is 

2 

expressed  in  m  . )  The  computational  expressions  used  for  G(z) 
and  the  dimensionless  quantity  Q(z)  are 


(B9) 


G(z)-9. 979(10"^)  +  9.382(10”®)a/^  +  9.900(10~^)z  , 


2 

in  (cmkHz)  with  z  in  m,  and 
(BIO)  [q(«)]-^  -  ti^o+n,o 

-  4.131  [i0"^)/gTzJ  +  (9.381)10~®//^jij[g(z)]"^  . 
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